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In this paper is shown an application of Clifford algebras to the construction of computationally 
universal sets of quantum gates for n-qubit systems. It is based on the well-known application of 
Lie algebras together with the especially simple commutation law for Clifford algebras, which states 
that all basic elements either commute or anticommute. 

PACS numbers: 03.67.Lx, 03.65.Fd, 02.20.Sv 



O 
O 

<N . 

Or 
<■ 
' 

(N ' 

> : 
■ 

o ■ 
o : 

o ■ 

o : 

^ ■ 

9 L,: 
^— > ■ 

c ■ 

s : 
cr 



X 



I. INTRODUCTION 



In this paper is discussed an algebraic approach to the 
construction of computationally universal sets of quan- 
tum gates. A quantum gate U for a system of n qubits is 
a unitary 2 n x 2™ matrix. It is possible to write U — e lH , 
where H is the Hermitian 2™ x 2™ matrix. 

A set of quantum gates Uk is (computationally) univer- 
sal if any unitary matrix can be obtained with given pre- 
cision as a product of matrices Uk- Algebraic conditions 
of universality can be described using the Lie algebra u 
of the Lie group of unitary matrices [QJ|]: if there is a 
set of Hermitian matrices and if it is possible to gen- 
erate a basis of space of all Hermitian 2™ x 2" matrices 
using only the commutators i[H, G] = i(HG-GH), then 
Uk = cxp(irHk) are a universal set of quantum gates if 
r is small enough. 

In this paper is presented an alternative approach to 
the construction of a universal set of gates using both Lie 
and Clifford algebras. It is possible because the algebra 
C(2™ x 2") of all 2" x 2" complex matrices is the com- 
plex Clifford algebra with 2n generators, i.e., there are 2n 
matrices Fk with the property /}} = /feP/ + i^Pfe = 
2Ski 1 (where 1 is unit matrix) and 2 2n different products 
of F fc generate a basis of C(2™ x 2") ||@. 

The 2n matrices I\ are not enough for proof of univer- 
sality, because we may not use arbitrary products of / & , 
but only commutators. In this paper it is shown that by 
using commutators of it is possible to generate only 
the (2n 2 + n)-dimensional subspace, but it is enough to 
add only one element r u and the new set is universal, 
i.e., it generates a full 4™-dimensional space u(2 n ). 

All 2n matrices may be chosen to be Hermitian 
and the full complex algebra was used for simplification. 
The extra Hermitian matrix is F u = 2/123 = ^1-^2-^3 or 
^1234) or any such product of three or four different TV 

A constructive proof of universality using the language 
of the Clifford algebras is based on a simple commuta- 
tion law of 4™ basic elements: they either commute or 
anticommute, because any such element is a product of 
up to 2n TV Direct construction of any 2" x 2™ ma- 
trix Jj = Jlfce/ °f the Clifford basis by commutators 
of 2n + 1 initial elements is shown below in Sec. II D, 
theorem 1. 



The question about universality is widely investigated 
|l|,||J|-P, but the method discussed in the present work 
has some special properties. Construction of a univer- 
sal set of gates uses only infinitesimal and continuous 
symmetries of group U(2") and does not require such 
discrete operations as permutations of qubits or basic 
vectors related to the "classical limit of quantum cir- 
cuits." The properties of discrete, binary transforma- 
tions of qubits simply emerge here from the structure 
of infinitesimal transformations of Hilbert space, i.e., di- 
rectly from Hamiltonians, cf. Elp|. 



II. CLIFFORD ALGEBRAS 

A. General definitions 

For n-dimensional vector space with a quadratic form 
(metric) g(x), the Clifford algebra 21 is a formal way to 
represent a square root of —g(x) ||,|| or, more formally, 
—g(x)t where 1 is the unit of algebra 21. The vector 
space corresponds to the n-dimensional subspace V of 

21: x 1— > x = J2?=o x i c h where x, e; € V C 21. From 
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x 2 = -g(x), i.e., (j2]Lo Wl) = T,i,jlo 9ijXi%j follow 
the main properties of the generators e; of the Clifford 
algebra: 



{ti,ej} = c,c, • c,c, = -2gij. 



(2.1) 



Let gij be diagonal and gu = ±1 (the case ga — is 
not considered here, but see H). Then, 



e,e 
e 
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= ±1. 



(2.2a) 
(2.2b) 



It is clearer from Eq. (2.2) that it is possible to gener- 
ate no more than 2™ different products of up to n C;. A 
linear span of all the products is a full algebra 21 [Q . Let 
us use the notations e/ = Ci 1 i 2 ...i k — ■ ■ • ti k , where k 
is the number of multipliers or the order of ej , k = A /"(J) . 

If there are no algebraic relations other than Eq. ( [2.2] ) , 
then the algebra has a maximal dimension 2 n and is 
called the universal Clifford algebra, £[(3), because for 
any other Clifford algebra 21 with the same metric g(x) 
there is a homomorphism £l(g) — > 21 (see Ref. 



1 



Le t us use the notation €1(1, m) for the diag onal m etric 
Eq. (2.2) with I pluses and m minuses in Eq. J2.2b ), i.e., 
for pseudo-Euclidean (Minkowski) space M. l ' m . There is 
a special notation for Euclidean space: €l(n) = €l(n, 0) 
and €[+(n) ee £[(0,n). 

Complexification of any Clifford algebra €1(1, rri) with 
I + m = n is the same complex algebra C), because 
all signs in Eq. ( 2.2b ) may be "adjusted" by the substi- 
tution e,t — » iefc. 

Let us denote ej ee y^ejej, i.e., if ef = 1, then ej = ej, 
but if ef = —1, then ef = iej and so always (ej) 2 = 1. 



B. Matrix representations 

All complex Clifford algebras in even dimension 
€l(2n, C) are isomorphic with a full algebra of 2™ x 2" 
complex matrices HQ- The simplest case €1(2, C) is the 
Pauli algebra. Matrices <7 X and <r y can be chosen as gen- 
erators eo, ei and a z is itoei — e^. 

The Pauli algebra is four-dimensional complex alge- 
bra and can also be considered as eight-dimensional real 
algebra, £[+(3). Prevalent applications of Clifford alge- 
bras in the theory of NMR quantum computation [ |U)|Jll| 
are based on a real representation £l + (3) rather than on 
a complex one £1(2, C), discussed in the present work. 
These two approaches are very close, but may be differ- 
ent in some of the details. 

There is simple recursive construction of the com- 
plex Clifford algebra with an even number of generators 
€l(2n, C) with €1(2, C). For n = 1, it is the Pauli algebra, 
and if there is some algebra €l(2n, C) for n > 1, then 



€l(2n + 2, C) €l(2n, C) <g) £1(2, C). 



(2.3) 



The proof of Eq. (2.3) is as follows: if eo, . . . , t2n-l are 
2n generators of £I(2n, C) then l 2n ® eo and l 2n ® ti 
together with 2n elements e^ (£> e^ are 2n + 2 generators 
of £[(2n + 2,C). 

Direct construction of €l(2n, C) is MM 



r 2k = t ® ■ ■ ■ ® t <%>a x ®a^ 

n—k— 1 

r 2 k+i = t<8>- -<g>t <S> (J y ® a z 

n—k— 1 



ut,, (2.4a) 
(2.4b) 



with e/ = 7], e 2 = 1, V/ € 0, . . . ,2n - 1. More gen- 
erally, algebraic properties of elements e; used in the 
paper are the same for different matrix representations 
e ; = MFiM- x , where M e SU(2"). 



C. Spin groups 

Most known physical applications of Clifford algebras 
are due to spin groups. The group has 2:f homomorphism 
with an orthogonal (or pseudo-orthogonal) group and is 



related to the Dirac equation B and the transformation 
properties of wave functions in quantum mechanics. 

E ach e lement x 6 V (see the definition of V above in 
Sec. II A) has an inverse x _1 = —x/g(x) if g(x) ^ 0. All 
possible products of even number of such elements with 
\g\ = 1 is the spin group. It is Spin(n) for €l(n) and for 
£[+(n). The group has 2:1 homomorphism with SO(n). 
For s £ Spin(n) an element of SO(n) is represented as 
r s : x h- > sxs -1 |f| . 

Because only products of an even number of elements 
of €l(n) are used in the definition of Spin(n), the group 
is a subset of even subalgebra €l e (n) c €i(n). In the Eu- 
clidean case, €l e (n) is isomorphic with €l(n — 1) and due 
to the property Spin(n + 1) may be defined as a subset 
of£l(n). 

Construction of the Spin(n + I) group from €l(n) 
is sometimes called the spoin group Spoin(n) = 
Spin(n + I). 

Let us consider (n + I)-dimensional space Al V, i.e., 
combinations y = A + x, x £ V. Let A(y) ee A 2 + g(x). 
The elements have an inverse (A + x) _1 = (A — x)/A(y) 
if A(y) ^ 0. Products of any number of such elements 
with |A| = I is the Spoin(n) group Q. 

The group Spoin(n) is 2:1 homomorphic with SO(n + 
I). For s £ Spoin(n), an element of SO(n + I) is repre- 
sented as r s :y i-> sy(s / y 1 , where y = y n + J2?=o Vi'i 
and (') is the algebra automorphism defined with basis 
elements as z'j = (-l^hj §. 



D. Lie algebras and Clifford algebras 

Clifford algebra is Lie algebra with respect to a bracket 
operation [a, b] = ab — ba 0]. Here we prove a result 
that is necessary for the construction of a universal set 
of gates. 

Theorem 1. Let €l(n,C) be the Clifford algebra and 
n be even. There are enough n generators k — 
0, . . . , n — I and any element e/ with Af(I) = 3 or 
Af(I) — 4 to generate elements of any order only using 
commutators of these n + 1 elements. 

A proof of this result has several steps. 

(i) If there are n elements eo, . . . , t n -\, it is possible 
by using commutators to generate also all elements of 



2z i z j 



2ti 



second order, i.e., [ej,ej] 

(ii) If there are all elements of second order and an ele- 
ment of third order, for example eoi2, it is possible to gen- 
erate any element of third order, i.e., 2eoim = [*oi2) *2m]j 

(iii) Analogously, if there is any element of order 2fc+l, 
it is possible to generate any element of the same order 
using no more than 2k + 1 commutators with elements 
tij ■ 

(iv) If we have all elements of third order, it is possible 
to generate any element of fourth order, 2tijki — [tijk, t{[. 

(v) Analogously, if we have all elements with the order 
M(I) = 2k + 1, it is possible to generate any element of 
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order 2k + 2, 2tj u i = [ej, ej], where I I. 

(vi) If we have an element of fourth order, it is possible 
to generate some element of third order, 2tijk = [tijkl , ti] 
(and so any element of third and fourth order). 

(vii) Analogously, if we have an element of order 2k+2, 
it is possible to generate some element of order 2k+l [and 
so any element with the order 2k + 1 or 2k + 2 as in the 
steps (hi) and (v)]. 

(viii) We have all elements with order less than or equal 
to 2k, k > 2 due to steps (i), (ii), and (iv) and we can 
prove the theorem by recursion: by using a commutator 
of an element with order 2k — 1 and an element with or- 
der 3 it is possible to generate an element of order 2k + 2 
and so any elements of order 2k + 1 or 2k + 2, as in the 
step (vii). 

Note 1. Instead of elements eo,...,e„_i, it is pos- 
sible to use eo together with n — 1 elements e; i : 

[eojCoi] = 2ei, . . . , [ej_i, = 2e/. 

Note 2. If n is odd, it is impossible to generate only an 
element with the order n, because due to step (vii) of re- 
cursion it would be generated only from an even element 
with the order n + 1, but there are no such elements. So 
in this case we need n + 2 elements, the extra one being 
eo,...,n-i- 

Note 3. If we use only n generators ti, then together 
with n(n — l)/2 commutators [efc,ty] = 2t)~j, k ^ j, it 
is possible to generate n + n{n — l)/2 = n(n + l)/2 el- 
ements, because, as may be checked directly, any new 
commutators may not generate an element with order 
more than 2. It is the Lie algebra of the Spoin(n) 
group, because products of exp(eefc) « 1 + eck belong to 
that group and the dimension of the group is the same, 
dimSpoin(n) =dimSO(n+l) = n(n + 1)/2. Despite the 
fact that only elements Ci, Af(I) < 2 belong to the Lie 
algebra, all 4™ elements Ci, Af(I) < n of <£l(n) belong to 
the Lie group Spoin(n) by definition and so a linear span 
of these elements is the full Clifford algebra. 

Note 4- The theorem was proved rather for the more 
general case of the Lie algebra of the complex Lie group 
GL(N, C), N = 2 n ' 2 of all matrices M, det(M) ^ 0, than 
for the unitary group U(N) C GL(iV,C). The proof for 
the Lie algebra u(N) of the unitary group U(iV) is di- 
rectly implied. It is sufficient to choose the initial ma- 
trices in u(7V) for a given representation, after which the 
Lie brackets may produce only matrices in u(N) for each 
step of the proof. 

It should be mentioned that there are two traditions for 
representations of u(iV). In physical applications, Hermi- 
tian matrices H are used, the Lie brackets are i[a, b], and 
the unitary matrices are represented as U — exp(— irH) 
due to relations with Hamiltonians and a quantum ver- 
sion of Poisson brackets [Q. In Eq. ( |2.4| ), elements 
ti = r u itox2, and e i23 (and iz kh see Note 1), i.e., all ef, 
are Hermitian. In more general mathematical applica- 
tions, u(iV) are skew-Hermitian matrices = —A and 
"i" multipliers are not present in the expressions for the 
commutators and the exponents B, because A = iH. 



III. APPLICATION TO QUANTUM GATES 

A. Universal set of quantum gates 

Now let us discuss the construction of universal gates 
more directly. Instead of Lie algebra u(2"), we should 
work with Lie group U(2"). Then an element corre- 
sponds to a unitary gate UJ = exp(iref ). One of the 
advantages of elements ej is the analytical expression for 
the exponent: 



UJ 



cos(r) + i sin(r)ej 



(3.1) 



Equation fl3.1| ) is valid for any operator with the property 
e 2 = 1 and it is true for all 4™ basis elements e?. 



It is also possible due to Eq. (3.1) to combine the ap- 
proach with infinitesimal parameters r jl],^| and an ap- 
proach with irrational parameters |5|||. The smaller r 
is, the higher is the precision in gen erat ion of arbitrary 
unitary gates in |Q,0|. Due to Eq. (3.1), accuracy may 



be arbitrarily high if we use gates Uj = e me ' with ir- 
rational tu/tt because for any r there exists the natural 
number N and e < r: Uf = (Ui) N . It should be men- 
tioned that the unitary gates do not necessarily have ir- 
rational coefficients even if w/tt is irrational, for example 
Uj = 0.8 + 0.6ef. 

Yet another advantage of the elements is a simpler 
expression for "commutator gate" . In the usual case jy,U , 
it is generated as 

e iri[H h ,Hi] _ e is/TH k ^s/TH, e -i^FH k & -iV7H t 



and the expression has precision 0(r 15 ). For elements 
ej, there is an exact construction. If Hi = t° and 
Hj = t°, then either [Hi, Hj] = or [Hi, Hj] = 2HiHj. 
The first case is trivial and for the second case due to 

Eq. (Eyh, 



e iri[Hi,Hj}/2 _ 6 -tHiHj _ ^Hi^tHj^i^H, 



After construction of the basis of Hermitian matrices 
Hi = e j , it is possible to use an expression 



/ \ N 



The expression has accuracy O (J2 a !/-^0- 

The approach to a universal set of gates U is more con- 
venient and constructive if we know the Hermitian matrix 
H, U T — e lrH . It is not a principal limitation, because 
for physical realizations we should know the Hamiltonian 
to construct the gates. It is also related to the universal 
quantum simulation Jj| in which H is the Hamiltonian 
and t is a real continuous parameter, the time of "appli- 
cation." 
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The description with an exponent may be even more 
complete, because by using H it is possible to find a 
unique U — exp(iif), but using U it is not always possi- 
ble to restore H because there are many H's for the same 
U . A simple example is U = ia a <8> 073 with two arbitrary 
Pauli matrices: U = (<r«®W®<r(>) = e'"" 8 "''. 



B. Two-qubit quantum gates 



Let us show how to build a universal set of one- and 



two-qubit gates using Eq. (2.4). For example, it may 
be 2n + 1 gates expire/), where cj are eo, with 
I = 1, . . . , 2n — 1, and ieoi2: 



r e 2fe,2fe+i 



;*2&+l,2fc+2 
1 



e = <g> a x , 

(g) C 

10( „_ fc _ 2) 



0~x 



eoi2 = 1 



5(n-2) 



(7., 



1 



(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 



with k = 0, . . . , n — 1 or n — 2. The elements were dis- 
cussed in Note 1, and it was shown that they generate 
the full Lie algebra u(2"). 



C. Nonuniversal set of quantum gates 



and (p,x) is the scalar product. The plane waves (j> v 
correspond to n particles with definite momenta. If 
O € SO (TV), then a transformation define d o n the ba- 
sis as So: <f) p — > <f)o P is a symmetry of Eq. ( |3.3|) . It is an 
analog of the classical transition between two configura- 
tions with the same total kinetic energy in "billiard balls" 
conservative logic . 

The general Dirac operator j| is the first-order differ- 
ential operator Dn — Y^hLq 1 ^ e k{9 / dxk) with a property 
D 2 N = — Ajy . If to use the Dirac operator for factoriza- 
tion of Eq. (ph, 



(3.4) 



(fD N - A)(Sjv + A)*(x , . . . , x N -i) = 0, 



then each component of \& is a solution of Eq. ( |3.3| ) 
and the action of the Spin(iV) group on \& corresponds 
to SO (TV) symmetry So described above and it has 
some analog in the classical physics of billiard balls. A 
Spoin(iV) group is represented less directly, but it can be 
considered as a symmetry between two stationary solu- 
tions with different total energies. 

The example above shows that it is possible to find 
some classical correspondence for elements e/, Af(I) = 2 
of the spin group and maybe for generators A/"(J) = 1 of 
the spoin group, but the special element with = 3 

does not have some allusion with classical physics. 



In an interesting question was raised, asking which 
sets of gates are not universal (and why). 

Products of gates U£ = e irefe = cos(r) + itk sin(r) gen- 
erate a group Spin(2n + 1) = Spoin(2n) c U(2 n ) due 
to Note 3. It is an interesting example of nonuniversal- 
ity when only one extra gate like e ITC ° 12 may produce a 
universal set with "an exponential improvement" from 
a subgroup dim Spoin(2rt) = n(2n -I- I) to a full group 
dim [7(2") = 2 2 ™. 

This result is more important if the extra gate e lTe > 
with Af(I) — 3 or Af(I) — 4 has a different physical na- 
ture from the gates with M{T) = I and M{T) = 2. It 
is not cl ear fr om Eq. ( |3.2| ) with the extra gate generated 
by Eq. ( 3. 2d ) is the simple onc-qubit gate. But this is 
not so for physical systems with natural Clifford and spin 
structure. 

A possible reason is the Schrodinger equation for 
n particles without interaction jL3|: ih(dip/dt) = 
J22=i(^a/nT>a)ip, or using m a = m and the Lapla- 
cian Ajv with N = un variables, it is possible to write 
for stationary solutions with total energy E, 



(Ajv + X 2 )i/j(x 0i • • ■ , xn-i) = 0, 



(3.3) 



where A = y / 2mE / h. Let the dimension of one particle 
motion be v = 2 for simplicity, N = 2n. 



Let us consider a full basis </> p (x) • 
space C of wave functions ip £ C 



i(p>x) 



on Hilbert 

N 
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